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ACOUSTIC SYSTEMS CONTAINING CURVED DUCT SECTIONS 

by W. Rostaf inski 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 
ABSTRACT 

The analysis of waves in bends in acoustical ducting of rec- 
tangular cross section is extended to the study of motion near 
discontinuities. This includes determination of the characteris- 
tics of the tangential and radial components of the non-propagating 
modes. It is established that attenuation of the non -propagating 
modes strongly depends on frequency and that, in general, the 
sharper the bend, the less attenuation may be expected. Evalu- 
ation of a bend’s impedance and of impedance -generated reflections 
is also presented in detail. 

INTRODUCTION 

Rectilinear acoustic systems consisting of sections of straight, 
hard-walled ducts of different diameters, or containing aperatures 
and similar constrictions, are characterized by propagation param- 
eters independent of frequency which can be described by relatively 
simple relations. Unfortunately, formulations derived for straight 
lines are inadequate for acoustic systems containing curved sections. 
The basic reason for this is that curved ducts constitute a dispersive 
medium; in other words, wave velocity in bends is a function of the 
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frequency. While in a straight duct only non-rigid walls create a 
dispersive medium for the progressing wave in its basic plane wave 
mode, in curved duct sections the zeroth mode will depend on fre- 
quency even with rigid walls and must be treated accordir^ly. The 
reason for this lies in the inability of curved ducts to carry plane 
waves. A plane wave at the inlet of the bend will become non-plane 
in the bend; for that matter, it begins to change to non-plane already 
in the straight duct feeding the bend. Consequently, in the bend, be- 
sides the progressing tangential component of particle velocity there 
will be more or less pronounces radial oscillations. The well-known 
illustration of the nature of plane waves advancing in tubes involves 
imagining an axial partition inside the tube. Insertion of such a par- 
tition in a straight duct will not modify the parameters of a plane wave. 
In a bent section, however, a partition concentric with the curved side 
walls will completely modify the original distribution of velocities and 
pressures. It will also affect the level of radial oscillations. The 
reason for this is that a partition forms hvo concentric bends less 
sharp than the original bend and creates essentially new boundary 
conditions. 

In previous investigations of the characteristics of waves moving in 
(1 2 ) 

curved bends ’ it was demonstrated that propagation in cylindrical 
bends, away from any discontinuity, is characterized by a lack of back- 
wards reflections resulting from curvature of the walls. Only onwards 

reflections are present. This very basic result of analytical work was 

( 3 ) 

recently verified ex]ieri mentally by Cummings. 


On the other hand, in 
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duct systems containing bends, the propagating waves will be accompa- 
nied by non- propagating modes originating at discontinuities. In Ref. 1 
the evanescent waves appearing at junctions between a straight duct 
and a cylindrical bend were studied in great detail (using an exact 
solution of the wave equation), but analysis was limited to long waves 
and very low frequencies. In a subsequent study (Ref. 2), higher fre- 
quencies and higher modes were considered, but in the steady state only; 
non -propagating evanescent waves requiring the use of imaginary - 
order Bessel functions were left out. The problem, however, is that 
the description of the motion of waves in duct systems cannot be com- 
plete without some understanding of the nature and extent of the non- 
propagating modes, which, as will be shown, may propagate quite 
far along the curved duct before dying out. It will also be shown that 
the non -propagating modes significantly modify the acoustic field at 
the bend’s inlet. Thus far. the basic mode in hard-walled ducting has 
been investigated quite extensively (Ref. 2). The results yield a good 
understanding of motion in typical industrial ducting. The higher 
modes have been given some attention but in view of the predominant 
role of the lowest mode and of its exclusivity in a relatively wide 
range of low frequencies, the higher modes seem to be of lesser 
importance in engineering. In most cases the higher modes, if at 
all excited, will belong to the evanescent, non -propagating category. 
Consequently, in Ref. 4 energy flow in curved ducts was analyzed in 
the basic mode only. The conclusion of this work was that an important 
attenuation of sound should be expected in sharp bends. 
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This paper will give the results of an analysis of the changes 
that occur as waves move in curved cylindrical ducting. The physi- 
cal system considered uses cylindrical coordinates (r, 0 ) and con- 
sists of a two dimensional circular bend with inner- and outer-wall 
radii Rj and R2 respectively. The walls are assun*ed perfectly 
rigid. The bend is connected to a straight duct forming a discontinu- 
ity at the junction. For study of the non-propagating modes a vibrat- 
ing piston will be substituted for the straight duct. The important 
subject of a bend’s acoustic impedance and impedance -generated 
reflections will be developed in some detail. Except for the solution 
of the characteristic equation for finding roots, the entire analysis 
is an exact solution of the wave equation. 

The Non-Propagating Modes 

In the study of the long waves that is for kR^ « 1, (Ref. 1), 
using a simplified series expression of the characteristic equation of 
wave motion in bends 

(akR )J’ (kRj) - J’ (akR.)J’ (kRJ . 0 (1) 

m m m ^ '^m 

where akR^ = kR2 it was possible to determine the sii^le real root 
yielding the propagating mode, as well as the infinite set of pure 
imaginary roots 

^ mTTi/ln a 1, 2, 3, . . . (2) 

corresponding to the non-propagating, evanescent modes. 

The evaluation of the motion of waves of higher frequencies 
(Ref. 2) was based on closed form solutions of the Bessel equation 
corresponding to orders - (m 4 1/2), m 1, 2, 3 
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This type of approach, however, does not allow for determination 
of the imaginary roots. Instead use must be made of special expres- 
sions for Bessell functions of imaginary order = ip_. Such 

■' m m 

(51 

series, first published in 1844 by Boole , were given in a more 
convenient form, and tabulated, by Buckens^®^ in 1963. Because of 
the little-known form of these series, they will be given here for 
convenience. 

Bessel functions of a purely imaginary order and of real argument 
are solutions of the equation 
+ (x^ + = 0 

The fundamental pair of solutions is defined by 
= A^(x) cos (fjL In x) - B^(x) sin (|i In x) 

= A^(x) sin (ji In x) + B^(x) cos ()i In x) 
where 





^ n=0 

and a and b are given by recurrence formulae 
’2n-2 


^2n 


^^(2n-2) 


4n(|i^ + n^) 



with 


"*^(2n-2) ~ *^^(2n-2) 

4n((j^ + n^) 

a = 1 and b^ = 0. 
o o 
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The most general solution of the wave equation in cylindrical coordinates 
(using kr = x) is 

4> = T/ 

m 

where, with Bessel functions of pure imaginary order but of real 
argument 

(0 (kr) = A G,,(kr) + B F,,(kr) 

^m m M m ju 

Differentiating </)^(kr) with respect to the argument along with a 
rigid wall boundary condition, gives 

G^(akRj) ■ F^(kRj) - G^(kRj) ’ F ’/akR^) 0 (3) 

This characteristic equation of motion of the non-propagating 
modes parallels equation 1. For very low frequencies, this character- 
istic equation yields the roots given by Eq. 2, derived in Ref. 1. For 
higher values of kR^, however, the roots do not fit into such a simple 
relation. These roots have been evaluated on a digital computer by 
programming Eq. 3. The program iteratively locates and converges 
to the roots (eigenvalues) of Eq. 3. The convergence criterion is 
accuracy to the eighth decimal place. 

The program is not long nor complicated. Because of its sim- 
plicity it can be written readily for any specific application. Here, for 
illustrative purposes only, two curved ducts characterized by a - R 2 / 

Rj - 2 and 4 have been analyzed and the results are given in Figs. 1 
and 2. For the duct with a 2, the roots can be approximated by 
Eq. 2 for kR^« 1. With higher frequencies the first root becomes very 
small, vanishing at kR^ ^ 77. The next root (m 2) vanishes at kR^ - 2n. 
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For a = 4 the root pattern is very similar but the small root zone is 
shifted to lower frequencies and all roots are smaller than in the case 
of a = 2. Approximate values of the roots for ducts with a = 4 can 
be calculated by Eq. 2 provided that kR^ < 0. 5. 

As a consequence of the appearance of small roots in the 
solution of the wave equation for non-propagating modes, the term 
exp(-|Lift' will remain unattenuated until G grows to large values such 
as 77 or greater. In other words there exist frequency ranges in 
which the non* propagating waves will be pronounced even well down- 
stream from discontinuities and the steady state distribution may not 
be reached in typical industrial bends of 45° or 90°. In such bends, 
the distribution of particle velocities and of pressure cannot be pre- 
dicted unless the evanescent waves are taken into consideration. 

Because with increasing ’’a” the general level of roots decreases, 
the sharper a bend the more the non-propagating modes will be sus- 
tained in their existence. Furtherm. re, since the evanescent waves 
may propagate far into the bend, any wall treatment designed to eliminate 
or reduce noise must be carefully designed to be effective; All depends 
on the degree of attenuation on non- propagating modes. At frequencies 
at which the non- propagating modes tend to propagate far into the bend, 
the higher propagating modes will appear if excited. 

To illustrate the point, propagation parameters have been cal- 
culated for a bend with radii a - 2 and for waves with frequencies 
corresponding to waves characterized by Bessel functions of order 
{m + 1/2) for m ranging from 0 to 6. The propagating modes of 
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these waves were described in Ref. 2. It may be worthwhile to 
restate here the values of the parameters of these waves for the 
lowest propagating mode: 

1.5 2.5 3.5 , 4.5 | 5.5 6.5, 

• I f i 

1.0115 1.663312. 2869 I 2 . 8817 3.4523 4.0065 
Tne calculated non -propagating, tangential and radial, vibrational 
velocities at two tangential positions: at 0 = 77/64, that is, very 
close to the surface of the piston generating waves at the bend's inlet, 
(and therefore in the region of discontinuity), and at ^ = 7r/4; are 
shown on Figs. 3 and 4. The vibrational velocities are nondimension- 
alized using v^, the vibrational velocity of the piston. For clarity 
in Fig. 3 the tangential waves are shown in three plots. In the range 
of low frequencies, for kR^ < 1.6633, the distribution of velocities 
near the piston is nearly that of a potential vortex but the attenuation is 
fast and by t/ 4 the vibrations become insignificant, reduced by two 
orders of magnitude. The case of kR^ = 2, 2869 is a particular case 
of very small vibrations throughout the bend, of unusual distribution 
and of very satisfactory attenuation. At frequencies corresponding 
to > 77, the attenuation is again very fast. The evanescent, 
radial vibrational velocities are shown in two plots - Fig. 4, For fre- 
quency parameters kR^ < 2. 2869 the radial oscillations are relatively 
small and attenuate fast, reaching negligibly small amplitudes by 77/4, 
At higher frequencies (except for kRj = 2. 8817) the radial component 
is even less pronounced, and it attenuates very fast. It will also be 
noticed that the radial evanescent waves may be either positive or neg- 
ative, depending on the frequency range. 


m 


L-R 


0.5 
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Figure 5 gives the tangential vibrational velocities v^ for three 
angular positions; at 71 /64, close to the piston, at 77/8, and at 
0 ^ «, Data for v^ = were taken from Ref. 2. The values of 
v^ at 7 t/ 64 and 77/8 were calculated by adding to (or subtracting 
from) (v^)^ values of the (v^)evanescent 

Fig. 3. Fxcept for the case of frequency parameter kR = 2. 8817, v^ 
acquires its final distribution relatively fast in duct cross sections. 

At 77/64 the distribution is quite pronounced near the two curved 
walls of the bend, but in the central part of the bend’s cross section 
v^ = 1, as it should be. However, at kR^ ^ 3. 4523, where the first 
root of the characteristic equation has already vanished, the distri* 
button of v^^ for all 0 seems to be very close to the final distri- 
bution at 0 - 

Acoustic pressure variations in the bend’s cross sections are 
shown on Fig. 6. The data are nondimensionalized using p^ the 
acoustic pressure at piston corresponding to vibrational velocity v^. 
The first two graphs, pertaining to very low frequencies, indicate a 
very uniform pressure distribution near the piston’s face. At all 
other frequencies the variation of pressure with the radius is pro- 
nounced. For kR^ - 3.4523 attenuation of the non- propagating 
modes takes place at f) ^ 77/64. 

The variation of pressure in the bends, as shown, cannot be 
calculated using information on tangential velocity distribution 
except for very low frequencies corresponding to kRj < 1, be- 
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cause the solution of the wave equation in evanescent modes is a 
series solution. Each term of the series pertains to its m’th root 
and has its own integration constant. Only when the first term of 
the expansion yields a satisfactory solution, as in the case of the very 
low frequencies, is such direct recalculation permitted. The calcu- 
lated distribui '• of pressure in bends matches well, in its general 

3 

geometry, the experimental data of Cummings . 

The Bend’s Acoustic Impedance and Transmission 
The specific acoustic impedance of duct bends 
z = pc(kRjA^)(r,/R^) 

has been derived in Ref. " This quantity, integrated and averaged 
over the duct’s cross section (of unit depth) gives 
= pc(kR^/2i^^)(a 4 1) = pc’ 

where the index c indicates that the expression refers to the cross 
section’s area, and where c’ is the phase velocity of the progressive 
wave at the bend’s centerline, as derived in Ref. 1. The dimension- 
less impedance for a duct of unll depth is defined by 
Zn - (kR^/2i/^^)(a + 1) - c’/c (4) 

The have been calculated for tlie lowest propagating mode and 
for several duct parameters a R2/HJ and for a range of f**e- 
quencies, and the results have been plotted on Fig. 7. Values of z^ 
for a relatively narrow duct (a - 1.5) are very close to unity, indicat- 
ing that the specific acoustic impedance of a duct of low aspect ratio 
is nearly the same as the specific acoustic impedance of a straight 
tube (z^ - pc) and consequently is independent of frequency. All 
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other characteristics on Fig. 7 become progressively steeper with 
increasing parameter a, that is, with increasing sharpness of bends. 

( However, there is always one particular value of the frequency pa- 

i 

ramater kR^ at which - pc in spite of the fact that the waves at 
these frequencies are far from being plane. In connection with this, 
it will be noted in Ref. 2 that at some frequency, distribution of 
tangential vibrational velocity shows a minimum of distortion at a 
frequency approximately corresponding to - pc. The consequence 
of z being generally different from pc will show itself in the sound- 
transmitting ability of bends. 

Assume waves advancing in two directions of n, the incident and 
the reflected wave trains. Expressions describing these waves, the 
acoustic pressure and the tangential velocity, are: 

p = -p 00/0t = +pikcfA^e“*^^ - A^e^^]F(kr)e^^ 

VjT = 1/r = +ii/|''Aje“‘^ + A^e^^ ]F(kr)/re^^ 

where F(kr) stands for the function of radial dependence of motion 
under consideration. A similar set of equations can be written for 
propa^iation of the transmitted wave in a straight duct connected to 
the bend. In this secono set of equations must be replaced by 
kx and the ex))ression for the axial velocity by / = ?<i>/?x. The 

A 

condition of c ontinuity of pressure and velocity at the junction between 
the bent (incident and reflected waves) and the straight (the trans- 
mitted wave) ducts at = x 0, (both assumed extending to infinity 
to simplify somewhat the model) is 
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^prop U nonprop 


prop 


I 
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nonprop 


0=0 



(5) 


where p and v include the incident and the reflected waves. Now 
evaluating the numerator and the denominator of the left hand side of 
Eq. 5, written in the form 


p (i + u , 'p + Po/P 
^prop i prop ^2 ^ 


/d + 
prop 


) 


it is verified tlu.t the variable terms in the parentheses decrease with 
frequency and become very small at low frequencies. Consequently 


Eq. 5 gives, ( ^\ 

j 


(kr/i»)(A. - A^)/(Aj + Aj.) = 1 

which, integrated and averaged over the cross section of the bend, 
yields 


(kRj/2i^)(a + 1) (A. - A^)/(A. A^) = 1 

or 


z„(Ai - A^)/(Aj 4 Ap - 1 

Now, the ratio of the amplitude of the reflected wave to the amplitude 
of the incident wave, that is, amplitude reflection coefficient r, is 
r A^/A. = (z^ - l)/(Zj^ +1) (6) 

2 

and the sound reflr^ction coefficient n? r , while the sound transrnis- 

r 

Sion coefficient will be ~ 1 - 

In the case of two bends (numbers 1 and 2) of different curvatures. 
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but of the same cross section, and connected in series: 

r = f z - z\//z + z\ 

(l " 2 ")/(l " 2 

The standard form of the standii^ wave radio SWR = (1 + r)/(l - r) 
as derived for straight ducts and pipes applies here without change. 
Physically, the situation is analogous to the acoustic phenomena at the 
boundary of two media. Equations 4 and 6 allow us to calculate reflec- 
tions at the boundary between a ’’dense medium” (bend) and a ’’less 
dense medium” (straight duct), or vice versa, depending on frequency. 
These reflections depend, of course, on the sharpness of the bend 
(parameter a). For the great majority of bends met in industrial 
applications, they will be quite low. Reflections appear because the 
impedance of curved ducts is different from the impedance of straight 
lines of the same cross section. For analytical purposes a straight 
duct- bend system can be replaced, as far as transmissivity of waves 
is concerned, by a single straight line, provided that the velocity of 
sound c' in the new straight duct is equal to c’ - c * (as if con- 
taiiiing some medium other than air). For very low frequencies and 
bends of arbitrary sharpness, the phase velocity wilt be calculated by 
methods of Ref, 1. For higher frequencies it will be convenient to 
use Fig. 7, Experimental information rn the transmissivity of bends 
may be found in Cummings' work. He gives data on the transmission 
coefficient of bends of a = 1. 59 and a - 10.2. The first bend tested 
in the frequency parameter range of 1. 06 < kR^ <5.3 had the trans- 
mission coefficient = 1 throughout the frequency range, which 
checks verv well against the analytical result (a^ - 1.0 - 4x10”'^), The 
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second bend, tested in the frequency parameter range 0. 068 < kRj < 
0. 34, had a transmissivity slightly less than 1, 0,975 on the average. 
This again checks well with the approximate theoretical r-^sult 
0.994) for kR, = 0.068 with z = 1.17. Furthermore, in the case 
of the second bend, at kR^ =0.2 Cummings obtained = 1. This 
could be anticipated because at this frequency parameter c = c’ and 
= 1. The selectively bumpy experimental curve may be explained 
by the high probability of frequency controlled evanescent waves in 
the zone of measurements and by a lack of uniformity in the trans- 
mission characteristic. 

The natural oscillations in ducts containir^ cylindrical bends may 
be best studied on systems open at one end and closed on the other. 

As a starting point, we may write that the input point impedance of a 
straight duct with a closed end (with the far end impedance infinite) 
and with negligible attenuation becomes a pure reactance 
X = -jpc cot (kl) 

Now, the well known equations given in some texts, for instance 
Ref. 7, for natural frequencies for a system of two straight ducts of 
lengths Ij and I 2 , of identical cross section, and connected in 
series, are: 

tan (klj) = cot (kl 2 ) for a closed far end 

tan (klj) = - tan (klg) for an open far end 

which will be satisfied if 

f(2n+l)7r/2 

kdj +ij) -{ 

V^Trn 


open duct 
closed duct 
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Now, if we replace duct number 2 by a bend, the equations for a 
closed end system will be tan (kl) = cot (p6) and 0 + kl = (2n + 1)7t/2 
where 0, in radians, refers to the angular position of the closed end 
of duct number 2. For very low frequencies, u- kR^ = ^^m ^'^^m 

where is the radius of the bend's centerline and s is the length of 
the bend along this center line. Consequently p9 = ks which reduces 
the expression for the bend to an expression for a straight line. For 
higher frequencies, the matter is not so simple. A more general ex- 
pression must be used which for a closed-end system is 
tan (klj) = Xctg(ykl^) where y = ^^^2 ^ some ratio of the 

acoustic specific impedances of the two elements in series in a 
ducting system. Expression tan(kl^) • tan (ykl^) = X, for y = 1, 
yields tan (kl^) = Vx; in general, however, a set of points must be found 
to satisfy this equation. We shall remember that for a narrow duct, up 
to and Including a = 1. 5, X ~ 1 and the solution is straightforward. 

CONCLUDING REMARKS 

Propagating and non-propagating modes of waves in the acoustic 
frequency range have been examined using a model, a system consist- 
ing of straight and carved duct sections. A method has been presented 
for determination of the non-propagating, evanescent modes generated 
at junctions, and of reflections due to impedance mismatch at junctions. 

It was indicated that each particular case of ducting must be examined 
separately. Only general formulations can be given for all shapes of 
ducting. Numerical values depend on the specific application and can- 
not be predicted without close analysis. 
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